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Abstract 

We discuss the theory of jet events in high-energy photon-proton interactions using 
a model which gives a good description of the data available on total inelastic ^p cross 
sections up to y/s=210 GeV. We show how to calculate the jet cross sections and 
jet multiplicities and give predictions for these quantities for energies appropriate for 
experiments at the HERA ep collider and for very high energy cosmic ray observations. 
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The hadronic structure of the photon and its interaction with hadrons are cur- 
rently subjects of considerable interest. The ZEUS and HI groups at the ep collider 
HERA have presented preliminary results on the 7p photoproduction cross section 
at 210 GeV center-of-mass energy which show the rise with energy typical of other 
hadronic interactions. It was shown by Gandhi and Sarcevic in a simple model 
that this rise can be used to discriminate between different sets of photon structure 
functions. The hadronic behavior of the photon at very high energies (~ 10 8 GeV 
photon energy) is also of interest for the theory of air showers triggered by cosmic 
ray photons |§. 

We have developed an improved QCD-based model of the photon interaction with 
nucleons f| which gives predictions for the total inelastic 7p cross section which agree 
well with the HERA data. The model has the interesting feature that the information 
needed to calculate a]£ el is largely determined by high energy nN scattering. In 
this paper, we will use the model to discuss the total jet cross section crjf t and the 
probabilities for multiple-jet events, and give our predictions for those quantities at 
HERA energies and the energies relevant for cosmic ray experiments. 

We begin by sketching our model for inelastic 7£> scattering. The theory is de- 
veloped in more detail in ||. The incoming physical photon state \j) p hys can be 
expressed as a superposition of the bare photon and the set of virtual hadronic states 
\m) which can be distinguished in the subsequent interaction, and appear with proba- 
bilities V m . Accordingly, the total inelastic 7p cross section is given in a semiclassical 
picture as 

°lnel = (1 - V had )O dlr +J2"Pm <T mp ', V had = ]T V m « 1. (1) 

m m 

Included among the states \m) are low- mass vector meson states such as the p, uj and 
0, and complex non-resonant final states which can be described on the average in a 
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quark-gluon basis. 

Whether the photon appears as a virtual vector meson state or a high-mass non- 
resonant state depends on the relative transverse momentum 2p±o of the quark and 
antiquark in the virtual transition 7 — > qq which initiates the hadronic interaction. 
We note that pj_ « 1/tj_, where r± is the average transverse separation of the virtual 
quarks during the lifetime of the qq system. If r± is greater than, or on the order of, 
the average transverse radius R± = l/Qo of a vector meson, QCD confinement effects 
will clearly set in, and the qq system will most likely appear in a hadronic collision 
as a light vector meson. For r_i_ < R± the qq system will be smaller than a vector 
meson. In this picture, the hadronic system behaves for p^ < Q like a vector meson, 
and for p±o > Qq, like a system of quasi-free quarks and gluons with a transverse 
area smaller than that of a vector meson by a factor (Qo/Rlo) 2 - Its interaction cross 
section will be reduced accordingly. 

After eikonalizing the hadronic cross sections using an impact parameter repre- 
sentation to account for possible multiple parton scatterings in a single 7p collision, 
the total inelastic 7p cross section is given by 0] as 

= ^r + W p Jd%(l-e- 2R ^ PP ) 

+ ^ e ;— / ^ U%{l-e- 2R ^ q ~ qP ). (2) 

q ft JQo P±Q J 

Here V p is the probability that a photon exists in the p meson state, V p = 47ta em / fl 2 , 
where j 2 p is the 7p coupling in the vector meson dominance model. We will treat 
the p, u and (f) as equivalent, and will use the quark model ratios for the couplings 
fv', then A = 4/3 for equal pp,ujp, and (ftp cross sections, and A = 10/9 for complete 
suppression of the </> contribution at low p±o ||. The third term in Eq. (2) includes the 
contributions from the excited hadronic states of the photon, which we collectively 
call qq states. The factor (e 2 q a em /'np 2 LQ )dp\Q is just the differential probability of 
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producing a qq pair at a relative transverse momentum 2p_i_o. 

The real part of the eikonal function for the pp interaction can be written in the 
form || 

Re X pp (b,s) = Re X Zft(s) + ReX%cD(b,s) 

= l^(b)[af oft (s) + a^ CD (s)}, (3) 

where A pp (b) is the density overlap function, 

A»(b) = J d 2 b' Pp (b) Pp (\b- b'\), J d 2 bA pp (b) = 1, (4) 

and (T^oft an d a QCD are the s °ft" an d hard-scattering parts of the intrinsic cross 
section. cr^ t was parametrized in f|] using a Regge-like form, 

<ft = + ^(a - mj)- 1 / 2 + a 2 (, - m 2 )- 1 , (5) 

while cr^f£, was identified with the inclusive parton-level cross section for pp scatter- 
ing, 

cr QCD = Yl 1 \ ,, / / rfx 2 ^ dp 2 ± ff(x 1 ,p 2 ± )f J p (x 2 ,p 2 ± )-^Y. (6) 

Using the equivalence of the p and 7r states in the quark model, we can equate the 
parton distribution and density profile functions for the p meson to the corresponding 
functions for the pion and take f[=f[ and p p (b)=p n (b). The parton distributions in 
the pion are known reasonably well [0j. The density profile functions p p (b) and p n (b) 
can be taken as the Fourier transforms of the electromagnetic form factors of the 
proton and the it meson. The density overlap function A pp « A 1 *" 9 is then given by || 

A pp (b) = ^^— 2 \vbK x (yb) - -^- 2 [K {ub) - K (pb)}\ , (7) 

47T p z — V A [ p Z — V A J 
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where K n (x) is a hyperbolic Bessel function, and the "size" parameters of the pion 
and the proton are p 2 = 0.47 GeV 2 and v 2 = 0.71 GeV 2 , respectively. The for- 
mer corresponds to a root-mean-square transverse radius for the pion or p meson 
R±=l/Qo « 2/ p. 

The eikonal function for the qqp scattering terms in Eq. (2) can be expressed 
similarly as 

2Re X ™ = A™(b, Pl _ )[a™} t (s,p xo ) + a^ D (s,p ±0 )}. (8) 

We assume that p q q{b) has the same form as p p (b) except with the size parameter 
fi = 2Q of the p meson replaced by 2p^ Q With this replacement, A qqp (b, p^ ) 
and pqq(b,p± ) are continuous with the corresponding functions for the p meson at 
Rlo = Qo, but describe a shrinking system for p± > Q . Finally, the cross sections 
a Toft ( s )Rlo) an d cr QCD(. s iP±p)y were assumed to scale with the physical cross sectional 
area of the qq system so that 

<fft(s,Pxo) + of&ohpu) = (Ql/plo)K P oft(s) + o& D (s)]. (9) 

This model guarantees that cr so / t and oqcd are continuous at p_i_o = Qo, and that the 
soft contributions from the qq states die out as l/pj_ as expected for higher twist 
contributions. 

We have calculated the total inelastic 7p cross section using Eq. (2) and the 
assumptions above. The soft cross section in Eq. (5) was determined || using the 
low-energy 'jp data 0. The calculation used the parton distributions of Owens [7J 
for the pion and those of Eichten et al. |K| for the proton. The parameter p± tm i n 



in Eq. (6) which determines the transverse momentum at which hard parton-level 
collisions come into play was taken from a similar fit to n^p scattering ||. The 
results agree with the preliminary HERA data |L| and those at lower energies as 
shown by the top curve in Fig. 1(a). The calculated rise in a]^ el at higher energies 



arises from the hard-scattering contributions, and is a prediction of the model rather 
than a fit to the HERA data. It is therefore of interest to look for other tests of the 
model. The jet cross sections predicted by the model provide one such test. 

The total jet cross section a]^ t (s,Q 2 ) is defined to be the part of the inelastic 
7P cross section which includes events with at least one semihard parton-parton (or 
7-parton) scattering with a momentum transfer p\ > Q 2 , irrespective of any soft 
processes that may occur. To find an expression for the jet cross section <7™f in 
the interaction between the proton and the hadronic state \m) of the photon, we 
use the fact that semiclassically exp[— 2RexQCD(b, s, Q 2 )] can be interpreted as the 
probability that there is no parton-parton scattering with p\ > Q 2 in a hadronic 
collision at impact parameter b. Using this observation, we can rewrite the expression 
for the total hadronic cross section a mp to separate out the jet-free part, 

= J d 2 b (l - e ~ 2Rex QCD - 2 Rex2 P f t^j 

= J d % ^ _ e -2Re X Z P ft '(b, s, Q 2 )^ e -2Re X % p CD (b, s, Q 2 ) 

+ J d% ( 1 _ e -2Re X ^ CD (b,s,Q 2 )y (10) 

Here 

a%( 8 , Q 2 ) = J d 2 b (l - e -2ReXQ P CD(b, s, Q 2 )^ (n) 

is the total cross section for events associated with \m) which contain jets with 
p± > Q. The factor <Jq CD which appears in the eikonal function in Eq. (11) is 
now to be evaluated using the expression in Eq. (6) with p\ min replaced by Q 2 . The 
remaining "nojet" cross section involves a modified eikonal function 

RexTfft'is, b, Q 2 ) = RexZt + R^XqU ~ R^dO, s, Q 2 ), (12) 
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and includes contributions from "soft" jets with p\ min < p\ < Q 2 as well as the 
usual soft term. 

The total jet cross section in jp scattering is the generalized sum of the cross 
sections cx™f (s, Q 2 ), weighted with the probabilities V m , over all the possible hadronic 
states of the photon, plus the very small jet contribution from <Tdir(s, Q 2 ), 

m 

= a dir (s,Q 2 )+XV p a^ et (s,Q 2 ) 

+Ee^L^<f(^ 2 ,pio)- (13) 

The jet cross sections for specific kinematic cuts, e.g., on jet angle, can be calculated 
using the same basic formulas with the cuts imposed on the integral for o"^, D in 
Eq. (6). 

We have calculated the total jet cross section a]^ t in the HERA energy range for Q 
= 2, 3, 4, and 5 GeV, using the same parameters as were used in ||. Our predictions 
are shown in Fig. 1 along with the calculated total inelastic 7p cross section and 
the experimental data ||. At ^=200 GeV, the jet cross section is predicted to 
be approximately 1.3% of the total inelastic cross section for Q= 5 GeV in rough 
agreement with earlier predictions ||, and 25% of the total for Q= 2 GeV. 

The rapid growth of the jet cross section with energy is associated with the growth 
of the parton distribution functions in the photon and the proton at the small values 
of x which become accessible at high energies. This growth also leads to an increasing 
probability of the scattering of more than one parton in a single 7p collision. For this 
reason, the inclusive parton- level 7p cross section (Tq' cd , which counts each parton 
collision separately, is not the same as <rjf t . The latter has the overcounting due to 
multiple independent scatterings removed in the eikonalization. We show the effect 
of the eikonalization up to HERA energies in Fig. 2. The effect is not very significant 
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for Q = 3, 4 and 5 GeV at y/s— 200 GeV, but causes about a 10% reduction in the 
jet cross section for Q = 2 GeV. 

In Fig. 2 we show crjf t for the cosmic ray energy range for the same choices of Q 
as in Fig. 1, together with the calculated cross sections aj p el and Oq CD for Q = 2 
GeV. At very high energies, the jet cross section clearly comprises a large part of the 
total inelastic cross section, e.g., at E 1 = 2 x 10 8 GeV (or y/s pa 2 x 10 4 GeV) crjf t 
constitutes 75% of ofrL for Q = 2 GeV. Even for Q = 5 GeV, crjf t gives nearly 
40% of alZei- Also, as more of the incident partons scatter, the effect of eikonalization 
becomes quite large as shown for Q— 2 GeV. 

To calculate the multiplicity of jets, we note that the average number of parton- 
parton scatterings with p± > Q in a single 7p collision at impact parameter b asso- 
ciated with the hadronic component \m) of the photon is 

n m (b, s, Q 2 ) = a™ p CD (s, Q 2 )A m ?(b) = 2Re X % p CD (b, s, Q 2 ). (14) 

Since the parton-parton scatterings in our model are independent, the probability of 
having n scatterings (2n jets) in such a hadronic collision has a Poisson distribution, 

P™ p (b,s,Q 2 ) = ^[n m (b,s,Q 2 )} n e- n m(b,s,Q 2 )^ (15) 

The probability P n (s,Q 2 ) of having n scatterings relative to all the inelastic events 
at an impact parameter b is then 

Pn(s,Q 2 ) = (^)-T^ 

m 

= (^ el r l j:^[d 2 bPr(b,s,Q 2 ) 

= K^^E^i / d 2 b[n m (b,s,Q 2 )] n e-^(b,s,Q 2 ) (16) 
for n > 1. Po(s, Q 2 ) is defined to be the probability of having a soft interaction with 
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no associated jets, 

P (s,Q 2 ) = (aZ^T.'Pm I d 2 b(l-e- 2Re ^oft\b,s,Q 2 )\ e -n m (b,s,Qi)^ 

m J \ J 

Note that we have neglected the very small two-jet contribution to these expressions 
from the direct QCD interaction of the photon with the partons in the proton. 

We have calculated P n for n = 0, 1, 2 and 3 at y/s— 200 GeV using Q — 2, 3, 
4, and 5 GeV. The results are shown in Table I. For Q = 5 GeV the probability of 
having a 2-jet event (a single parton scattering) is about 1%; the probability for 2n-jet 
events with n larger than 1 is essentially zero. For Q = 2 GeV the total probability of 
having any hadronic jet event is approximately 27%, with the 2-jet events comprising 
more than 90 % of the total. 

We also investigated jet production in 7p scattering at very high energies. In 
Fig. 3 we show the distribution of P n at y/s = 2 x 10 4 GeV with the choice of Q 
= 3 GeV. It is more likely at this energy to have jets in an event than not, with 
an approximately 70% jet probability. While the 2-jet production is still the most 
probable, events with more than 2 jets occur with about 20% probability. 

It will be quite interesting to see if the predictions above are verified in future 
experiments since they connect the hadronic properties of the photon quite directly 
to those of the pion. In particular, the jet structure of pion and photon-induced 
reactions on protons should be quite similar. 
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FIGURE CAPTIONS 

FIG. 1. (a) The calculated total jet cross sections as functions of ^/s and the 
transverse momentum cutoff Q. The total 7p inelastic cross section calculated in B 
is compared to the low-energy || and HERA [Q data in the upper curve, (b) The 
eikonalized and inclusive jet cross sections for the same values of Q. 

FIG. 2. The total inelastic and jet cross sections at cosmic ray energies. o~qcd * s 
the inclusive jet cross section. 

FIG. 3. The probability distributions for n parton-parton collisions or 2n jets in 
a 7£> collision. 
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Q (GeV) P Pi P2 P 3 

2 0.734 0.228 0.020 0.002 

3 0.904 0.078 0.002 0.000 

4 0.954 0.030 0.000 0.000 

5 0.971 0.013 0.000 0.000 



Table 1: Probabilities P n of having n parton-parton collisions or 2n jets in a single 
inelastic 7p event at ^=200 GeV. 
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